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Abstract 

We study theoretically electronic transport through a contact of a quantum wire with 2D or 
3D leads and find that if the contact is not smooth and adiabatic then the conduction is strongly 
suppressed below a threshold voltage Vr, while above Vr the dc current I is accompanied by 
coherent oscillations of frequency / = J/e. The effect is related to interelectronic repulsion and 
interaction of dc current with the Friedel oscillations near a sharp contact. In short conducting 
channels of length L < Lq ~ hvp/eVr and at high temperatures T > To ~ eVr/kB the effect is 
destroyed by fluctuations. 

PACS numbers: 73.63.-b, 73.23.-b, 73.63.Fg, 72.10.Fk 
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In contrast with 2D and 3D systems where basic electronic properties are usually well 
described in terms of the Fermi liquid and single-electron noninteracting quasipariticles, ID 
systems of interacting electrons are better described in terms of the Luttinger liquid (LL) 
with bosonic excitations. The LL is an alternative to the Fermi liquid in ID (for a review 
see Ref. 1). Interaction in ID systems greatly affects both the electronic structure and 
transport. In particular, a power-law suppression of density of states arises near the Fermi 
energy and even isolated impurities strongly suppress conduction resulting in power-law 
dependence of conductance on voltage and/or temperature 1 . This behavior was confirmed 
experimentally in various ID systems including semiconductor quantum wires 2 and carbon 
nanotubes 3 , and was described in terms of macroscopic tunneling between different minima 
of a periodic potential tilted by external bias. The periodic potential is associated with the 
Friedel oscillations (FO) induced by an impurity. It was shown recently 4 that the power-law 
regime takes place only at small enough voltages, while above a threshold voltage a dynamical 
regime of conduction sets in. In this regime the dc current is accompanied by oscillations 
of frequency / = I/e. The effect is induced by motion of the FO in a repulsive electronic 
system when dc current passes an impurity. Interesting effects can also be expected at 
sharp, non-adiabatic, contacts between quantum wires and electrodes of higher dimension, 
i.e. between the Fermi liquid and strongly correlated electronic state. This problem has 
not been studied yet, and usually the boundary conditions derived by Egger and Grabert 5 
for ideal adiabatic contacts are used. These conditions were derived only for expectation 
values and, therefore, cannot be used to describe fluctuations, which are very important in 
ID systems. Further, real contacts are not necessarily adiabatic and one can expect that 
reflections of electrons from contacts may result in the FO and, hence, in effects similar to 
those predicted for the case of impurities. Here we derive boundary conditions for formation 
of the FO near contacts and show that the conductance is affected by the FOs, resulting 
in the dynamic regime of conduction that resembles the Josephson effect and the Coulomb 
blockade. 

Below we set e, h and ks to unity, restoring dimensional units in final expressions when 
necessary. 

We derive boundary conditions using the ideas of the scattering approach (for a review 
see Ref. 6) We describe a ID conductor as a potential barrier at \x\ < 1/2 with a channel 
forming a quantum wire along the rr-axis. The wire is attached to two symmetric 2D or 3D 
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metallic leads. The region of conducting wire is considered as a scatterer. Since the results 
are very similar for both contacts we consider here, for brevity, the left lead only and the 
result for the right one will be given without derivation. 

Without the loss of generality we assume that longitudinal (along the x-axis) and trans- 
verse motions are separable. The longitudinal motion in the leads is characterized by wave 
vectors k and energy E\ — J^, and transverse motion is described by energy e n , the total 
energy being e = E\ + e n , where n is an index labeling transverse energies. We assume that 
electrons in the leads do not interact. Then we solve an equation of motion for electronic 
field operators in the leads using the continuity of both the field operators and their deriva- 
tives at |x| = 1/2. This allows us to express the solution for the n-th transverse eigenstate 
in terms of the field operator tp b at the boundary 

%b(x) = tbb cos kx + -rd x ijjb sin kx. (1) 
k 

This expression contains both incident and outgoing waves. According to the causality 
principle, the incident wave ipi n (x) is determined by a state of the lead far away from the 
barrier. Therefore, tpi n (x) must not depend on properties of the barrier. Equating the 
incoming part of Eq. (1) to the form describing free particles we find 

*Pb - jrd x ip b = -^L J2 d n,k5(e - e n - ^—), (2) 



where ki = ^2m(e — e n ) and c n ^ is an annihilation operator of an electron in the lead with 
a longitudinal momentum k in the n-th transverse mode of the lead. 

Eq. (2) relates the field operator at the boundary to the equilibrium states of the n-th 
transverse mode of energy e in the lead. We are interested in finding a relation between the 
boundary value of the field operator corresponding to the lowest transverse eigenstate of the 
conducting wire and the incident state of the lead. To find this relation, we project Eq. (2) 
onto the eigenstates of the wire. Since transverse states of the lead are not eigenstates of 
the wire, we obtain an infinite system of linear equations for boundary values of the field 
operators tpj of the transverse eigenstates j of the wire 

'•'./ -'Y.'' J r i) (3) 

r 

We must find a solution of Eqs. (3) for the state j = describing the lowest subband which 
is responsible for an electronic transport in the wire, while the states j > with higher 



transverse energies do not contribute to the transport. It follows from Eq. (3) that the 
relation we are looking for has a form 

A(e)^ + B(e)dJ = ]T l(k)c n 2n5(e - e n ), (4) 

where the exact expressions for the coefficients in Eq. (4) depend on the shape of the contacts. 
The boundary condition for the right contact has the same form but with complex-conjugate 
coefficients. 

Coefficients in Eq. (4) are not arbitrary. In particular, they must provide correct anticom- 
mutation relations for electronic field operators. It is worthwhile to relate the coefficients to 
such physical parameters of the system as transmission probability t of incident electrons. 
Therefore, we consider the system of noninteracting electrons, for which we can easily solve 
the equations for the field operators inside the wire. Then we impose a requirement of 
fulfillment of anticommutation relations, and calculate the conductance. This allows us to 
reduce the number of undetermined constants. As it is more convenient to express bound- 
ary conditions in terms of physical values, we multiply Eq. (4) on the left by its Hermitian 
conjugate, then we transform the obtained equation to the time representation assuming 
that the coefficients are slowly varying functions of energy in the region close to the Fermi 
energy. Finally, we find boundary condition for the left (right) contact 

V -fp ± j + v F fp F = 1 £ ctc^-** (5) 

n,n' 

where j and p are operators of current and of the smooth part of charge density perturba- 
tions, pp is the 2/c^-component of charge density, which is related to the FO, / is a number 
ofthe order unity if the transmission probability is not close to unity, and / ~ ^2(1 — t) if 
1 — t <C I. Thus the FOs disappear if the contacts are adiabatic. 

In order to check the validity of conditions (5), we considered a wire with noninteracting 
ID electrons with smoothly widening contacts, so that the contacts are nearly adiabatic. 
We also assumed that there might be a potential step of the height Uq <C Sf at the interface. 
Under these assumptions we were able to use the quasiclassical approximation in the lead 
and match the quasiclassical solution outside the ID conductor with the exact solution inside 
the channel. We found that the condition (5) yielded a correct result for the conductance 
G = tG in agreement with the Landauer formula. Here G = e 2 /h is the conductance 
quantum. 
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In order to take into account interaction in the quantum wire, we consider spinless (spin- 
polarized) electrons described by a bosonic displacement field $ p obeying the Tomonaga- 
Luttinger Hamiltonian 1 . The bosonic field $ p determines current and perturbations of charge 
density by means of relations 

1 <9<f> L 

p = — + — cos (2k F x - 2$) , / = (e/n)dt$. (6) 

it ox ix 

The interaction is assumed to be short-range, described by the parameter K p < 1 character- 
izing the strength of interaction (K p = 1 for noninteracting electrons). A short-range interac- 
tion corresponds to gated quantum wires where the long-range part of interaction is screened 



by 3D gate electrodes. For quantum wires we can roughly estimate K p ~ yhv F e/e 2 ~ 0.2 
^/et>ir(cm/s)/10 7 , where e is the background dielectric constant. 
In the Heisenberg representation $ p satisfies the wave equation 1 

(v 2 dl-d?)%(t,x) = (7) 

where v = v F /K p is the velocity of plasma waves. This equation must be solved with 
boundary conditions at the contact which we obtain after bosonization of (5). We should 
note that since one assumes a linear dispersion of electrons within the LL theory, the theory is 
valid provided that all energies are small in comparison with the Fermi energy. However, the 
term which is responsible for the FO after bosonization has the form ^ e F cos(2l> p +k F l). 
Generally, this value is of order of the Fermi energy and it is small only if \J\ — t <C 1. 
Therefore, we consider nearly adiabatic contacts, where y/1 — t <C 1. Further, following 
Ref. 5 we take into account screening of the potential of the leads by a 3D gate. Finally, we 
obtain the boundary conditions for bosonic field $p at the left and right contacts in the form 

^d x $ p =F d t $ p + fe F cos(2$ p =F k F l) = P L , R . (8) 

Now we represent the bosonic field operator as a sum of its expectation value and a fluc- 
tuating part, = § p {x = =F//2) = + <£l,r- Then we perform thermodynamic 
averaging of both sides of Eq. (8) and obtain boundary conditions for the expectation values 

j^d x <p L)R =F d t <p LyR + d sin(20 Lji? + a) = U LtR , (9) 
p 

where U F , R is a potential applied to the left (right) contact, d and a are given by expres- 
sions d = fe F \J (cos 2ip) 2 + (sin 2ip) 2 , a = arctan fe^gf) f° r eacn con tact. The equation for 



the expectation values 4>l,r which determines current and smooth perturbations of charge 
density contains the term related to the FO. This term depends on fluctuations in the ID 
channel. If the contact is adiabatic / = this term disappears, and Eq. (9) is reduced to 
boundary conditions by Egger and Grabert 5 . 

In order to find correlation functions for the bosonic fields, we need to find first the cor- 
relation functions for fluctuating parts of the operators 8Pl,r = Pl,r — (Pl,r) ■ In frequency 
representation it reads 

{{6P L (u),6P L (u/)}) =4ir 2 ucoth^5(tu + iu'). (10) 

Now we can calculate the current induced by voltage V = Ur — Ul applied to the leads. If the 
contacts are adiabatic then / = ^2(1 — t) = 0, and there is no FO at the contacts. In this 
case, we obtain an ohmic current j = d t &/n = G V. The result is different if the contacts 
are not adiabatic. In this case, we solve Eq. (7) by performing the Fourier transformation 
and substituting then the solution into the boundary conditions given by Eq. (8). In this 
way we derive the equations for the field operators &l,r at the corresponding contacts. 

A(cj)$ l (w) + B(u)$ R (w) + S L (w) = P L (u), 

B(u)$ l (uj) + A(u)$ R (u) + S R (u) = P R (u), (11) 

where A(u) = u(i - KpS l nuJti ), B(u) = ucotut h t x = l/v, S l , r (uj) = 
fe F J dte luJt cos(2$ iifi (t) =f k F l). We cannot solve these nonlinear equations easily, the main 
difficulty being the account of fluctuations. We assume that fluctuations are Gaussian. 
Strictly speaking, the fluctuations are not Gaussian. However, our approach can be justified 
strictly in the case of strong interelectronic repulsion and in the limit of high voltages, where 
the non-Gaussian part of fluctuations is small. This can be shown similarly to the case of a 
current passing through an impurity in a ID conductor 7 . 

We consider several limiting cases which can be solved analytically. First, we try to find 
a stationary solution in the case of low applied potentials ±V/2. After averaging Eqs. (11), 
we obtain the following equation for each contact 

dsin(20 Lji? + a) = V/2. (12) 

Eq. (12) has stationary solutions for a finite voltage when d ^ 0. We can calculate d by using 
self-consistent harmonic approximation 1 , in which fluctuations are assumed to be Gaussian. 
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In this approximation, we replace sin2y3 with e 2 ^2y3, and obtain a simple expression 

d = 2fe F e~ 2 ^. (13) 

We have obtained linear equations for fluctuations which can be solved easily. Thus we can 
find <£l,r and using the expression for anticommutators of fluctuation sources (10) we can 
calculate the mean square of fluctuations 

In pure Luttinger liquid without FOs at the contacts, this integral diverges logarithmically 
both at high and low frequencies. At high frequencies it must be cut off at the energies of 
the order of the Fermi energy. The divergence of fluctuations at low frequencies is a feature 
of one-dimensional systems. In our case we obtain that the integral is cut off at the lower 
limit of the order of d. As the complete expression for (<PL tR (u))) is rather cumbersome we 
give the result only for frequencies oj > d , which determine a large logarithmic contribution 
to (<Pl,r)- Since the result is the same for both the contacts, we omit the indices L, R 

7rcoth^ir p 2 [2-(l-K p 2 )sin 2 ^] 

{ip {uj)) = u 4iq + (i-ag)Wurtr (15) 

In addition to the logarithmically divergent part, this expression contains the oscillating 
factor induced by reflections of fluctuations from contacts. If the length of the quantum 
wire is large enough, / 3> v/d, these oscillations contribute little to the integral and the 
oscillating factor can be replaced with its average value K p /(1 + K p ). Further, the result 
of integration depends on a relation between d and temperature T. At low temperatures 
T <^ d, the integration yields with a logarithmic accuracy 

/^ = ^P_ ]n _Ei_ m (i 6 ) 

xy 1 1 + K P rfcos20 v ; 

Since d depends on (y? 2 ), Eq. (16) is self-consistency condition for (y3 2 ). Substituting d 
from Eq. (13), we find d and the maximal value of the left hand side of Eq. (12). Thus we 
determine the value of the threshold voltage for which a static solution for mean phase 
exists. 

V T ~ f—e F . (17) 

We see that in the case of interelectronic repulsion, K p < 1, the mean square of fluctuations 
(ip 2 ) and the amplitude of the FO are finite, while in the case of noninteracting system, when 



Kp — 1, fluctuations become infinite and d — 0. Note that the role of the FOs here is similar 
to that of impurities in LL. Thus we obtain that if the voltage applied to the contacts is 
small enough, V < V T) current does not flow through contacts. This result is a consequence 
of the approximation in which only Gaussian fluctuations were taken into account. If we 
took into account fluctuations of solitonic type for which the phase increases by 2it due to 
tunneling, we would obtain a small tunneling current at V < Vt- Fluctuations of this type 
were studied in the case of impurities and they resulted in power-law I-V curves 1 . 

Now, we consider the case T > d. The self-consistency equation has solutions which 



correspond to a finite value of fluctuations only if T < T ~ f 1 - K pe F , so at temperatures 
T > T Q the FOs are destroyed by thermal fluctuations and do not affect electronic transport. 

In the case of a short enough channel I <C v/d we must not average Eq. (15) over 
oscillations at outi < 1. At these frequencies (0 2 (u)) in Eq. (15) is proportional to u^ 1 as 
before but with a different factor. As a consequence the FOs do not affect the conduction 
when I <C v/Vt- 

If we increase the voltage above Vr, a transition to a nonstationary regime of conduction 
starts. We consider the limiting cases of low temperatures T <^ Vr and of long wire I ^> 
v/Vt- It is difficult to obtain I-V curves at low voltages accurately, the main difficulty 
being the account of fluctuations with the mean square value periodically depending on 
time. The problem is simplified at high voltages V 3> Vt when the mean square value 
becomes nearly constant with small oscillating component. In this case Eqs. (11) can be 
solved perturbatively assuming that the oscillating part is small both in the fluctuations 
(0 2 ) and in the mean phase (0). Note that Eqs. (11) describe two interacting nonlinear 
oscillators. In such systems different solutions can exist. We consider a solution for which 
oscillations at both contacts are synchronous. Then the averaged Eqs. (11) can be reduced 
to a single equation, which for the constant applied voltage reads 



To solve this equation we need to calculate the value of d(t) which is determined by fluc- 
tuations. In order to do this we solve Eqs. (11) for fluctuations substituting d(t) in the 
form d(t) = do + d c cosuot + d s smu>ot, and assume that d c ,d s <C d. Thus we obtain a 
system of linear equations for <Pl,r, from which linear equations for correlation functions 
{{ip RtL (u), <Pr,l{u')}) can be easily derived. After a rather cumbersome but not too difficult 




(18) 
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derivation we find that the main logarithmic contribution to d is determined by Eq. (15) as 
before, but with different infrared cut-off frequency b ~ d 2 /cj <C d. Thus the self-consistency 
equation is again of the form (16), but with b instead of <icos20. From the self-consistency 
condition we find 

d = V T [ — , d c ~ — , d s ~ — In— . (19) 

We see that at high voltages the solution with finite amplitude of the FOs exists only when 
K p < 1/3, i.e. when interelectronic interaction is strong enough. The result differs from the 
stationary case where fluctuations do not destroy FOs at any repulsion strength K p < 1. 
The result also differs from the case of impurity, where the critical value at high voltages is 
K p = 1/2. 

Now we can solve Eq. (18) easily in the limit of high voltages V 3> V T) and to calculate 
current using Eq. (6). The total current consists of dc part, I = VG — I n i, and ac part, 
J ac sinco> t, which oscillates with frequency oj q = 2nl/e ~ eV/h 

iitv>i , + „„2 uoU T _ n _ (V T 



I ac ~ 27rG d K 2 p /(K 2 p + tan 2 -^), I nl ~ 2vr J I ac . (20) 

The oscillating factor in these expressions is due to reflections of generated current pulses 
from the contacts. 
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